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Abstract 



Materials such as foams, concentrated emulsions, dense suspensions or colloidal gels, are yield stress fluids. Their steady flow 
behavior, characterized by standard rheometric techniques, is usually modeled by a Herschel-Bulkley law. The emergence of 
Q-i techniques that allow the measurement of their local flow properties (velocity and volume fraction fields) has led to observe new 
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complex behaviors. It was shown that many of these materials exhibit shear banding in a homogeneous shear stress field, which 



cannot be accounted for by the standard steady-state constitutive laws of simple yield stress fluids. In some cases, it was also 
observed that the velocity fields under various conditions cannot be modeled with a single constitutive law and that nonlocal 
models are needed to describe the flows. Doubt may then be cast on any macroscopic characterization of such systems, and one 
i— —i may wonder if any material behaves in some conditions as a Herschel-Bulkley material. In this paper, we address the question of 
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the existence of a simple yield stress fluid behavior. We first review experimental results from the literature and we point out the 
main factors (physical properties, experimental procedure) at the origin of flow inhomogeneities and nonlocal effects. It leads us 
to propose a well-defined procedure to ensure that steady-state bulk properties of the materials are studied. We use this procedure 
to investigate yield stress fluid flows with MRI techniques. We focus on nonthixotropic dense suspensions of soft particles (foams, 
concentrated emulsions, Carbopol gels). We show that, as long as they are studied in a wide (as compared to the size of the material 
mesoscopic elements) gap geometry, these materials behave as 'simple yield stress fluids': they are homogeneous, they do not 
exhibit steady-state shear banding, and their steady flow behavior in simple shear can be modeled by a local continuous monotonic 
constitutive equation which accounts for flows in various conditions and matches the macroscopic response. 

Keywords: Yield stress fluid; Viscoplastic flow; Local rheology; Herschel-Bulkley 
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1. Introduction 

Materials such as dense suspensions, colloidal gels, microgel 
suspensions, concentrated emulsions or foams, are yield stress 
fluids: below a yield stress T y , they do not flowQ. Their steady 
flow behavior in simple shear is characterized in standard The- 
ological experiments, e.g. in a Couette or a cone-and-plate 
geometry QJJ] . It is usually well accounted for by a Herschel- 
Bulkley equation t(%) = T y + J] HB j", the parameters of which 
depend on the details of the material microstructure [2]. The 
ability of such models to help predicting flows in more com- 
plex configurations, e.g. in extrusion flows [4], then relies on 
the robustness of the characterization in simple shear. 

However, the recent emergence of techniques J2, Etl that al- 
low the measurement of local flow properties (velocity and vol- 
ume fraction fields) has led to revisit the behavior of pasty ma- 
terials. The reason is that, in rheometric experiments, one only 
measures macroscopic quantities (torque T(D.) vs. rotational 
velocity Q). The material constitutive law relating local quan- 
tities (shear stress r(y) vs. shear rate y) is then derived under 
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1 We neglect slow creep flows driven by coarsening dynamics in foams and 
emulsions Q]. 



the assumption that the flow and the material are homogeneous, 
and that the material can be modeled as a continuum with a 
local constitutive law. Local measurements have shown that 
these conditions are actually not met for some materials. The 
observations of steady flow inhomogeneities (shear banding) in 
homogeneous stress fields have in particular shown that, in ad- 
dition to their yield stress r v , some materials are also charac- 
terized by a critical shear rate y c below which they cannot flow 
steadily QG|, LZk. LSD - Moreover, in some materials, it was found that 
flows under various conditions cannot be described by a single 
local constitutive law J^, [Icjl, which led to propose nonlocal 
modeling of their behavior [9, . These observations contrast 
with standard laws inferred from purely macroscopic measure- 
ments, such as that of Herschel-Bulkley, which are local and 
allow steady flow at any low shear rate j. 

Doubt may finally a priori be cast on any macroscopic char- 
acterization of yield stress fluids, and one may wonder if any 
material behaves in some conditions as a Herschel-Bulkley ma- 
terial. The question we address in this paper is whether such a 
simple yield stress fluid behavior exists or not. In Sec. |2l we 
first focus on the shear banding issue; then, in Sec. [3] we in- 
vestigate the consistency between yield stress fluids flows and 
a local modeling of their behavior. Both sections are structured 
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in the same way. We first review experimental studies of yield 
stress fluid flows in the literature: we present the main sources 
of shear banding (Sec. 12. U and of discrepancy between macro- 
scopic and local behavior (Sec. 13. U . This leads us to distin- 
guish three classes of materials: thixotropic materials (such as 
colloidal gels), nonthixotropic dense suspensions of rigid par- 
ticles (such as noncolloidal suspensions and colloidal glasses), 
and nonthixotropic dense suspensions of soft particles (such as 
foams, concentrated emulsions and Carbopol gels). We discuss 
the role of thixotropy and shear-induced migration in the stud- 
ied complex behaviors for the first two classes of materials. We 
then focus on nonthixotropic dense suspensions of soft particles 
(Sees. 12.21 & 13. U , whose case is yet unclear as contradictory 
observations are reported in the literature. We first show that 
the above mechanisms should not be at play in these materials. 
We then point out issues related to the procedure used to study 
their flows. We propose a well-defined procedure to ensure that 
steady-state bulk properties of the materials are studied, and we 
investigate the flows of several foams, concentrated emulsions, 
and a Carbopol gel, with MRI techniques (Sees.l231&ET2l We 
show that, as long as they are studied in a wide (as compared 
to the size of the material mesoscopic elements) gap geometry, 
these materials behave as 'simple yield stress fluids': they are 
homogeneous, they do not exhibit steady-state shear banding, 
and their steady flow behavior in simple shear can be modeled 
by a local continuous monotonic constitutive equation which 
accounts for their flows in various conditions and matches their 
macroscopic response. 

2. Shear banding 

In this section, we address the question of shear banding in 
yield stress fluids. We discuss the phenomenon and the ori- 
gins identified in two classes of such materials: structural inho- 
mogeneities in thixotropic systems; volume fraction inhomo- 
geneities in dense suspensions of rigid particles. We then in- 
vestigate flow stability in nonthixotropic dense suspensions of 
soft particles. We discuss the relevance of the above mecha- 
nisms of shear banding in these systems. We use a procedure 
designed to circumvent possible artefacts and show that these 
materials do not display any steady-state shear banding, in con- 
trast with several results of the literature. Possible explanations 
of these discrepancies are proposed at the end of the section. 

2.1. Phenomenon and possible origins 

Shear banding in yield stress fluid flows here refers to the 
coexistence of liquid-like (sheared) and solid-like (unsheared) 
bands at steady-state in a homogeneous shear stress field! It 
is observed in some materials when a macroscopic shear rate 
Tmacro lower than a critical value y c , which depends on the mate- 
rial, is applied. In these conditions, the material basically splits 
into a region flowing at % and a non flowing region, whose rel- 
ative extent ensures that the shear rate spatial average equals 



■/macro 1 12 , 0] • The critical shear rate y c is the shear rate 
below which no stable flow is possible. It appears to be a new 
mechanical characteristics of these systems, in addition to their 
yield stress r,,: at steady state, flow implies both r > r v and 
y > j c - In an inhomogeneous stress field, flow localization oc- 
curs with any yield stress fluid; the signature of shear banding 
is then that the material flows at a nonzero shear rate y c at the 
interface between the flowing region (where the shear stress t 
is higher than t v ) and the non-flowing region (where t < t v ) 
J3l. This phenomenon is generally attributed to the fact that the 
theoretical material steady-state flow curve is non-monotonic 

Q. 

Although shear-banding materials cannot flow steadily at a 
local shear rate smaller than y c , any macroscopic shear rate 
tmacro < % can be applied to such system by the relative motion 
of two boundaries. This implies that their macroscopic charac- 
terization may fail to represent their actual behavior. Indeed, 
typical constitutive laws T(yrn aC ro) found near the yield stress 
in macroscopic experiments are such that y macro tends continu- 
ously to zero as t approaches t v ; this is the case of standard laws 
such as the Herschel-Bulkley law. This contrasts with the lo- 
cally observed behavior of shear-banding materials where yi oca i 
tends towards a critical shear rate y c + as r approaches T y , 
This issue is illustrated in Fig.Q] 




"It is also sometimes called 'discontinuous' shear banding to avoid confu- 
sion with other kinds of flow inhomogeneities. 
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Figure 1 : Illustration of the difference between the material intrinsic constitu- 
tive law (grey line) and the flow curve measured macroscopically in a homoge- 
neous stress field (black line) in a shear-banding material. At high macroscopic 
shear rate (y ma cro > 7c)> fl° w is homogeneous and the macroscopic flow curve 
matches the local behavior. At low macroscopic shear rate (y macl0 < y c ), shear 
banding occurs: the local shear rate in the flowing material —%- differs from 
the macroscopic shear rate and the macroscopic flow curve does not match the 
local behavior. 



Identifying which materials exhibit shear bands and which do 
not is thus of high importance; only the last ones may behave 
like simple yield stress fluids and be characterized using stan- 
dard rheometric experiments. Two possible origins of shear- 
banding in yield stress fluids have already been identified exper- 
imentally: structural inhomogeneities in thixotropic systems, 
and volume fraction inhomogeneities in dense suspensions of 
rigid particles. 
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2.1.1. Structural inhomogeneities in thixotropic systems 

Shear banding in yield stress fluids seems to be observed 
in all thixotropic colloidal gels 1 12 , 15 , 13 , KJ 01 ■ In these 



materials, it is attributed to the competition between aggrega- 
tion, which is due to attractive interactions and thermally ac- 
tivated structuration, and shear, which tends to break the ag- 
gregates. Starting from a liquid state, when the applied shear 
rate y is too low, the material basically splits between a liquid- 
like region (where it is a suspension of independent aggre- 
gates) and a region that gets continuously structured under 
shear until it jams (i.e., a percolated network forms). Shear 
banding is also observed in dense suspensions of rigid noncol- 
loidal particles, where it is attributed to competition between 
sedimentation-induced contact formation and shear-induced re- 
suspension jl7lll8ll . When the local shear rate is too low, the 
shear resuspension mechanism is not sufficiently efficient to 
prevent contact formation between the particles and the mate- 
rial jams as a percolated network of particles form. 

In these two kinds of systems, shear banding seems to result 
from shear/structure coupling, and from the development of two 
bands of different structures, similarly to what is observed in 
wormlike micellar systems B19I1 . Several models based on the 
idea that the material properties depend on its structure and that 
the structure depends on shear history lEol 21 , 22 ] have indeed 



predicted shear banding in simple shear, as a result of the insta- 
bility of some solutions of homogeneous structure, and describe 
correctly the flow inhomogeneities observed in thixotropic sys- 
tems |2]. Such models can be an explanation for shear banding 
in a given material only if a macroscopic thixotropic behavior - 
i.e., dependent on shear history i23ll - is found experimentally. 



2.1.2. Volume fraction inhomogeneities in suspensions of rigid 
particles 

Colloidal glasses are not much thixotropic: they get struc- 
tured at rest and destructured at flow start-up, but their flow be- 
havior does not depend significantly on their flow history (see 
e.g. pil). The above mechanism is thus not expected to be at 
play in such systems. However, shear banding has recently been 
observed in these materials and explained to result from the de- 
velopment of an inhomogeneous volume fraction profile [25]. 
The basic idea is that a small change in the volume fraction may 
result in a significant change in the rheological properties, in 
particular in the yield stress value, which is an increasing func- 
tion of <p. In a homogeneous stress field, at low applied shear 
rate, regions of lower <p (and thus of lower r y ) may then flow 
while regions of higher <p (and thus of higher r y ) are jammed. 
The development of an inhomogeneous volume fraction profile 
in a homogeneous stress field suggests an instability induced 
by a shear-concentration coupling [26, 25]. Banding resulting 
from volume fraction inhomogeneities has also been observed 
in dense rigid noncolloidal suspensions [fl7l|27ll, where shear- 
induced migration 11281 l29l 13011 creates jammed regions of vol- 
ume fraction <p higher than the jamming packing fraction <p m . 

It is worth noting that these systems are made up of rigid par- 
ticles. Indeed, although hydrodynamic interactions may play 
a role, recent theoretical developments lf30ll supported by ex- 
perimental findings [27] show that the shear-induced migration 



efficiency is strongly enhanced by direct rigid contact forces be- 
tween the particles. This characteristic may thus play a crucial 
role in the emergence of shear banding in these systems. 

2.2. Nonthixotropic dense suspensions of soft particles 

The case of nonthixotropic dense suspensions of soft parti- 
cles (foams, concentrated emulsions, microgels) is less clear, as 
contradictory observations are reported in the literature. Becu et 
al. 113 ill reported shear banding in adhesive emulsions, whereas 
Ovarlez et al. I32I1 did not observe any shear banding in the 
same system as well as in four other adhesive and nonadhe- 



sive emulsions. Rodts et al. 113 311 re ported shear bands in 3D 
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sheared foams, but Ovarlez et al. 113411 did not observe any shear 
banding in the same system and in three other foams. In 2D 
foams (bubble rafts) Gilbreth et al. I35T1 reported shear banding 
whereas Katgert et al. II 1 CW1 did not. This suggests that subtle 
differences in the materials, setup, or procedures, may exist, 
and have to be identified. In the following, we first discuss the 
possible relevance of the mechanisms mentioned in Sec. l2.ll in 
these systems. 

2.2.1. Nonthixotropic behavior 

The systems we deal with are basically nonthixotropic: their 
flow behavior does not show significant history dependence. 
This is illustrated in Fig. [2^, where the response of an emul- 
sion, a foam, and a Carbopol gel to consecutive up and down 
stress ramps is displayed. Such experiment is classically used 
to assess the thixotropic behavior of a material [2]. During 
the up-ramp, as long as the shear stress is below the (static) 
yield stress, the material is strained in its solid regime; this 
explains why the strain rate is roughly constant and propor- 
tional to the ramp rate in a first stage fl. Once the material 
flows, there is no significant difference between the responses to 
the up- and down-ramps: these materials can be considered as 
nonthixotropic. This contrasts with the hysteretic response ob- 
served in a thixotropic material (bentonite suspension, Fig. I2J?), 
which is the signature of strong dependence of the behavior on 
flow history, i.e., of thixotropy. 

This suggests that the mechanisms discussed in Sec. 12.1.11 
should not be at play in the Carbopol gel, the foam and the 
emulsion, as their structure should remain homogeneous when 
they flow. Structuration at rest is nevertheless observed in some 
of these systems ll39l lii PloTl . i.e., their structure in their solid 
and liquid regimes may differ. E.g., after a preshear, foams pro- 
gressively recover their elasticity at rest (G' increases in time) 
as coarsening induces structural rearrangements i39l. I3I1 . Struc- 
turation at rest can lead to differences between the static and 
dynamic yield stress 1401 . As detailed in Sec. 12.3.11 this can 
lead to transient flow inhomogeneities but not to steady-state 
bands. 

Recently, models of soft-jammed systems such as concen- 
trated emulsions (a nonlocal model 12211 and the SGR model 
1 21]) have been extended to account for possible bulk shear 
banding. Both models predict macroscopic shear banding asso- 
ciated with structural inhomogeneities. We stress that, in these 
models, the behavior depends on internal variables that depend 
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Figure 2: Shear stress t vs. shear rate y when consecutive up (symbols) and down (lines) stress ramps are applied, with a constant rate a, (a) to a foam (up triangles 
and dashed line, a=\ Pa/s), a Carbopol gel (empty circles and full line, a=0.5 Pa/s) and a concentrated emulsion (empty diamonds and dotted line, a=0. 1 Pa/s), 
and (b) to a thixotropic (bentonite) suspension (a=0.5 Pa/s). A 4°angle cone-and-plate geometry is used, with serrated surfaces (60 mm diameter) for the foam, 
with sandblasted surfaces (40 mm diameter) for the other materials. The emulsion is composed of 1 /im diameter dodecane droplets dispersed at a 73% volume 
fraction in water stabilized by Sodium Dodecyl Sulfate at a 1 wt% concentration (see [36] for details). The Carbopol gel is made up of Carbopol 980 dispersed and 
neutralized at a 0.3% volume fraction in water (see |37] for details). The foam (Gillette Regular) consists of 27 /jm diameter bubbles at a 92% volume fraction. The 
bentonite suspension is made up of (smectite) clay particles of length of order 1 /jm and thickness 10 nm suspended at a 6% volume fraction in water 13811 . 



on shear history, i.e., the system is supposed to be thixotropic: 
they cannot apply to nonthixotropic suspensions of soft parti- 
cles. 

2.2.2. Material homogeneity 

As discussed in Sec. 12.1.21 volume fraction inhomogeneities 
may lead to shear banding. One may wonder whether they can 
be responsible for some of the reported shear bands in dense 
suspensions of soft particles. Shear-induced migration is well 
documented for rigid particles ll28l I29I 17 , 30[ 23] . Migration 



of deformable particles, and particularly drops, has been much 
less studied. It is found that single droplets migrate away from 
rigid walls, due to asymmetric flows around deformed droplets 
1I42II43I1 . but binary collisions tend to homogenize the systems 
for0 > 10%. 

To test the possible emergence of volume fraction inhomo- 
geneities in dense suspensions of soft particles, we study the 
evolution in time of the bubble or droplet volume fraction pro- 
files when foams or concentrated emulsions are sheared in an 
inhomogeneous stress field (where particles migrate towards 
the low shear zones in suspensions of hard spheres). We use 
a wide gap Couette geometry, where the shear stress distribu- 
tion is: 



t{R) = T(RdR*/R 2 



(1) 



where R is the radial position in the gap and R/ is the inner 
cylinder radius. The materials and methods are described in 
detail in Ovarlez et al. J^Sl. The average bubble diameter of 
the studied foam is 73 /mi, its gas volume fraction is 89% [34]. 
Five emulsions are studied; their droplet diameter and droplet 
volume fraction are: (0.3 jum,75%), (1 //m,75%), (6.5 yum,75%) 
and (40 ,um,88%); both adhesive and nonadhesive emulsions 
are studied H32I1 . The inner cylinder radius of the Couette cell 
is 4.1 cm; the gap is 1.9 cm wide. The systems are sheared 
at various constant rotational velocities Q, ranging from 1 to 



100 rpm, during long times (from 1 to 24 h, depending on the 
system). The Couette cell is inserted in a Magnetic Resonance 
Imager 1I44I1 and droplet or bubble volume fraction profiles are 
measured as in [ 17, 32]. 
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Figure 3: Concentration profiles observed within the gap of a Couette geometry 
by MRI techniques after shearing a concentrated adhesive emulsion (squares, 
average applied strain: 70000), and a foam (open circles, average applied strain: 
12000). The emulsion is composed of 6.5 /im diameter silicone oil droplets 
dispersed at a 75% volume fraction in a mixture of 50 wt% glycerine and 50 
wt% water stabilized by Brij and trimethyl tetradecyl ammonium bromide at a 
6.5 wt% concentration 1 32]. The foam is composed of 73 /im diameter bubbles 
at a 89% volume fraction in a SLES foaming solution described in |34]. 



Typical results are shown in Fig. [3] We find that concentrated 
emulsions and foams remain homogeneous within +0.2%, even 
when large strains are applied; in the experiments shown in 
Fig. [3] strains are larger than 10000. The same result is found 
with all emulsions studied in 113211 : no other foam is studied. 
This is in contrast to the behavior of dense suspensions of rigid 
particles. In particular, the strain to reach a stationary inhomo- 
geneous state can be as low as 20 in the suspension of rigid 
particles studied in [27], of particle size similar to the droplet 
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and bubble size of an emulsion and the foam we have studied. 
The absence of observable migration in concentrated emulsions 
and foams may thus be due to the very different nature of inter- 
particle interactions in these systems: indeed, Lhuillier lf30h has 
recently pointed out the major role of direct rigid contact forces 
in migration. Other factors may play a role; in particular, we 
emphasize that the dispersed phase volume fraction of all the 
soft materials studied here is well above the jamming transi- 
tion, which may hinder local increase of the particle volume 
fraction. 

From the absence of observable migration for very large 
strains, we conclude that this process, if it exists, would be 
much too slow and inefficient to lead to observable shear band- 
ing in dense suspensions of soft particles at the usual experi- 
ment timescales. 

Finally, for nonthixotropic dense suspensions of soft parti- 
cles, thixotropy and migration should not be at the origin of 
shear banding. To explain the apparently contradictory obser- 
vations of the literature, one may thus wonder whether a yet un- 
known mechanism is at play in some of these systems only, or 
if differences in the setup or procedure are involved. In the fol- 
lowing, we present an experimental investigation of flow stabil- 
ity in various materials with a well-defined procedure designed 
to circumvent possible artefacts. 

2.3. Investigation of flow stability in dense suspensions of soft 
particles 

We propose to study the possible shear banding in non- 
thixotropic dense suspensions of soft particles as quantified by 
an intrinsic bulk property, the critical shear rate below which 
flow is unstable. In order to design appropriate setup and pro- 
cedure, we first deal with possible artefacts that may lead to ob- 
serve flow inhomogeneities that are not steady-state bulk shear 
bands. 

2.3.1. Possible artefacts 

Finite size effects. In systems of small gap to dispersed ele- 
ments size ratio, structural inhomogeneities are sometimes ob- 
served. Indeed, the local shear stress / shear rate relationship 
r(y) in concentrated emulsions 45] and 2D foams ifioll in 
confined geometries is found to depend on the system size and 
on the boundary conditions, suggesting that nonlocal laws are 
needed to describe their behavior J^l [Tl|] (this is discussed in 
Sec. 0. The material's apparent viscosity r/y near the walls 
for a given applied stress can in particular be smaller than that 
of the bulk material |0,|45t]: the material is locally 'fluidized' 
by a rough surface. This mechanism could a priori lead to ob- 
serve coexisting liquid-like and solid-like bands (note however 
that no shear-rate discontinuity is reported in |0, 1(1 45]); this 
would not reflect a bulk property. 

Transient flow inhomogeneities. Divoux et al. J46ll have re- 
cently observed long-lived transient flow inhomogeneities dur- 
ing the flows of Carbopol gels in a thin gap Couette cell. These 
inhomogeneities always disappear at steady state l46ll : at any 
rate, steady-state flow is homogeneous. However, given the 



very long lifetime of these inhomogeneities at low shear rate 
(they may last as long as 1 day at a 0.1 s rate), one may won- 
der if all bands reported in the literature are really steady-state 
bands. It is thus crucial to understand in which condition such 
transient inhomogeneities are generated; this may allow us to 
propose a procedure to circumvent this possible artefact. 
It is worth noting that a slight structuration at rest is observed 
in the studied systems |46[ 4(J: their elastic modulus increases 
with the resting time, and their static and dynamic yield stress 
differ. Moreover, in each experiment, flow is imposed on the 
material initially at rest: the transient bands characterize flow 
start-up. A first basic explanation for the development of these 
bands can be given in the framework of a simple model, illus- 
trated in Fig.Hk. Let us consider a material characterized by a 
simple yield stress fluid behavior r = f(y) in its liquid regime, 
with f(y) — > Ty d when y — > 0, which defines its dynamic yield 
stress T yd . In addition, we suppose that flow can be induced on 
the material initially at rest only when t > T ys , with r ys > T yd , 
which defines its static yield stress r Vi ; this model accounts for 
the possibility of a specific structuration at rest. If a macro- 
scopic shear rate y maC ro is applied on the material initially at 
rest in a slightly inhomogeneous shear stress field (as in a thin 
gap Couette geometry), two cases are found, depending on the 
value of y ma cro- If /(Tmacro) ^ T y , after the solid/liquid tran- 
sition, the material simply flows homogeneously at the applied 
shear rate y maC ro and the resulting stress is r = /(y m acro)- If 
/(tmacro) < Ty s (Fig. H^), although the material should be able 
to flow at a stress lower than the static yield stress, the stress 
applied on the material in its solid regime has first to reach T yi . 
Consequently, just after the solid/liquid transition, shear must 
localize in a band in the region of higher stress (near the in- 
ner cylinder in a Couette cell): the local shear rate in the band 
tband nas to t> e higher than the macroscopic shear rate y maC ro 
to ensure that the resulting stress r = /(yband) equals t Vi , i.e., 
tband = / 'll^)- Note that the mechanism proposed by Moor- 
croft et al. 14711 in the framework of the SGR model is basically 
based on the same idea. 

Of course, it is only a rough explanation of the Divoux et al. 
|46] observation^. However, the above simplified model prob- 
ably captures the basic reason for the existence of a banded 
flow when flow is initiated on a structured state: stresses at flow 
start-up have to be higher than in steady state. It also suggests 
that other loading procedures can be used to avoid generating 
transient bands; in the following, we propose such procedure 
(Fig.Hp). 

2.3.2. Setup and procedure 

We are here interested in steady-state shear bands. As de- 
fined in Sec. 12.11 these should result from flow instability at 
low applied rate. We stress that, to investigate their possible 
development, a homogeneous flow of the material in its liq- 
uid regime should first be imposed at high shear rate; next, the 



It does not explain what happens subsequently to the band, i.e., how the 
material at rest is progressively eroded by the flowing band. Moreover, the full 
story of flow initiation is far more complex, as shown by Divoux et al. 14011 ; in 
particular, flow is initiated by full slippage at the rotor of the Couette cell just 
after the solid/liquid transition. 
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shear rate should be decreased continuously (see Fig. \4fc). In 
these conditions, transient flow inhomogeneities due to specific 
structuration at rest (Sec. I2.3.11 l should not be observed: this 
procedure is thus preferable to applying a given shear rate on 
the material initially at rest, which may lead to wrongly report 
shear bands if experiments are not long enough (we recall that 
transient bands can last for 1 day). Then, if the material is not a 
shear-banding material, homogeneous flow should be observed 
at any low rate. In contrast, if the material is a shear-banding 
material, flow instability should develop at low shear rate. Note 
that this procedure also avoids observing the history-dependent 
flow inhomogeneities recently predicted by Cheddadi etal. [ 48 ] 
for materials exhibiting significant normal stress differences; 
this point is discussed in more detail below. 

We focus on the behavior of dense suspensions of soft par- 
ticles: concentrated emulsions, Carbopol gel, and foams, de- 
scribed in detail in J32 , 4§l 34 ] . To better show the difference 



between these materials and thixotropic systems, we also study 
the behavior of a colloidal gel (bentonite suspension), described 
in J3l ■ The materials are sheared in the same wide gap Couette 
geometry as in Sec. 12.2.21 This ensures that the bulk behav- 
ior, free from the finite size effects discussed in Sec. 12.3.11 is 
studied (this is shown in detail in Sec.[3]l. Consistently with the 
above procedure, we first apply a high rotational velocity Q. to 
the inner cylinder of the Couette cell, and we decrease progres- 
sively the value of Q. Steady-state azimuthal velocity profiles 
V(R) are then measured as a function of the rad ial p osition R in 
the gap using MRI techniques described in 1501 14411 . for various 
constant rotational velocities £1 

2.3.3. Experimental results and analysis 

All the studied materials exhibit similar behavior; an exam- 
ple is shown in Fig. [5] We observe that the whole sample is 
sheared at high rotational velocity Q. When decreasing Q, be- 



low a critical value that depends on the material and is of order 
40 rpm in Fig. [5] the material is sheared only in a fraction of 
the gap: V(R) vanishes (within the measurement uncertainty) 
at some radius R C (Q.) inside the gap, which decreases when Q. 
decreases. 




4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 

Radius R (cm) 

Figure 5: Dimensionless velocity profiles V(R)/V(Rj), measured by MRI tech- 
niques, for the steady flows of a concentrated emulsion in a wide gap Couette 
geometry, at various rotational velocities ranging from 5 to 130 rpm (see leg- 
end); the solid line is the theoretical profile for a Newtonian fluid. The emulsion 
is composed of 1 /im diameter silicone oil droplets dispersed at a 75% volume 
fraction in water stabilized by Sodium Dodecyl Sulfate at a 8.5 wt% concentra- 
tion 132 1. Figure from Ovarlez et at 1331 . 



Flow inhomogeneity here originates from the shear stress 
field inhomogeneity (Eq.Q]). Shear localization should then not 
be here confused with the shear banding observed in homoge- 
neous stress fields: it is a feature of any yield stress fluid flow 
in a Couette geometry. At low D., flow stops at a radius R c 
within the gap where the local shear stress t(R) equals r y (i.e., 
R c = Ri Jr(Ri)/Ty). The decrease of R C (Q.) as Q is decreased 
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then arises from the rate dependence of the constitutive law at 
the approach of t v . 

Although it induces shear localization for both shear-banding 
and non-shear-banding materials, the stress inhomogeneity can 
be used advantageously to analyze the possible shear-banding 
behavior from the shape of the localized velocity profiles. In- 
deed, from the constitutive behavior point of view, the ques- 
tion of shear banding can be posed as follow (Sec. 12. U : what 
is the value of the shear rate as t tends towards the dynamic 
yield stress? In shear-banding materials, it should tend towards 
a nonzero value y c ; in simple yield stress fluids, it should tend 
towards zero. When flow is localized, r approaches t v when R 
tends to R r . As the local shear rate is 



y(R) = VI R - dV/dfl 



(2) 



the local velocity in a shear-banding material should thus tend 
to zero with a nonzero slope |dV/dR| = y c at the interface be- 
tween the sheared and the unsheared regions, as in a homo- 
geneous stress field U7[]. In contrast, the slope of the velocity 
profile at the interface between the sheared and the unsheared 
regions should be equal to zero in a non-shear-banding material. 

In the above analysis, we have assumed that the yield cri- 
terion only involves the shear stress t, whose value thus de- 
termines the position of the interface between the sheared and 
the unsheared regions. If the material displays significant nor- 
mal stress differences in its solid regime (as e.g. in dry foams 
lf5lln . those may enter in the yield criterion and things can be 
different. Whereas steady-state normal stresses in the flowing 
region depend only on the shear rate and are thus independent 
of the initial conditions, different nonzero normal stresses can 
be stored in the material at rest in the nonflowing region: their 
value depends on shear history, i.e. on preparation. Some 
sample histories may then lead to a discontinuity of the nor- 
mal stress differences at the interface between the sheared and 
unsheared regions [48]. When taken into account in the yield 
criterion, this discontinuity implies that the shear stress at this 
interface is higher than in the case of continuous normal stress 
differences and that y at this interface is nonzero B48I1 . This 
leads to a slope discontinuity in the velocity profile that appears 
as a critical shear rate which is not an intrinsic property of the 
material, but rather reflects the loading procedure. We point out 
that such discontinuities are not expected to take place with the 
loading procedure we use, as the solid-like zone is progressively 
formed by decreasing continuously the inner cylinder rotational 
velocity. We thus think that even with materials exhibiting sig- 
nificant normal stress differences, only intrinsic shear banding 
will be evidenced. 

To analyze the possible shear banding, we now plot localized 
velocity profiles obtained for various materials, and character- 
ized by almost the same position R c (here ^5.4 cm) of the inter- 
face between the sheared and the unsheared regions (Fig. |6^). 
We focus on the shape of the velocity profiles at this interface. 
From this plot, it clearly appears that (i) the thixotropic sys- 
tem (bentonite suspension) exhibits an abrupt transition from 
flow to rest, as expected from the shear-banding behavior of 
these materials observed in cone-and-plate geometry lfl2i @], 



whereas (ii) concentrated emulsion, foam and Carbopol gel ex- 
hibit a smooth transition from flow to rest and thus seem to be 
not shear-banding materials. 

To better illustrate these differences and their link with the 
constitutive behavior of the materials, we combine the local 
stress data t(R) (Eq.[T} and the local shear rate data y(R) (Eq.|2]l 
extracted from the localized velocity profiles (Fig. |6^,b); this 
allows the reconstruction of the local constitutive law r(y) at 
the proximity of the yield stress r y . These laws are depicted in 
Fig. [(J. A clear difference between the thixotropic system and 
the nonthixotropic dense suspensions of soft particles is again 
observed: when t — > t v , y tends towards a finite value - the 
critical shear rate y c (here of order 6 s _1 ) - in the bentonite 
suspension, and towards zero in dense suspensions of soft par- 
ticles. We have observed this smooth transition from flow to 
rest, with the same procedure, in several other emulsions [32] 
and foams 113411 . 

We conclude that none of the studied dense suspensions of 
soft particles is a shear-banding material. Of course, due to fi- 
nite experimental resolution, we can only provide upper bounds 
for the value of the critical shear rate - if any. Note that the ac- 
curacy of the shear rate measurement depends on the spatial 
resolution and on the MRI signal to noise ratio, which depends 
itself on the material and on external factors. We finally ob- 
tain typical upper bound ranging between 0.01 s and 0.3 s _1 



2.4. Probable explanation of previous observations 

From MRI observations on various systems, we claim that 
nonthixotropic dense suspensions of soft particles do not show 
any bulk steady-state shear banding. To strengthen this claim, 
we have to check that previous observations of shear rate dis- 
continuities in similar systems in the literature can be explained. 
In the following, we review these observations and identify pos- 
sible origins of the discrepancy between these results and ours. 



Finite size effects. Denkov et al. 115 211 have reported shear band- 
ing in 3D foams, in a parallel plate geometry of gap to bubble 
size ratio of order 10. Gilbreth et al. 113511 have reported shear 
banding in bubble rafts, in a geometry of gap to bubble size 
ratio of order 20. Such ratios are probably much too small to 
ensure that a bulk property is studied (see Sec. |3). It is worth 
noting that Katgert et al. 11 1 CHI reported nonlocal effect in bubble 
rafts of similar gap to bubble size ratio as in [35]. It is thus pos- 
sible that the II351 15211 shear rate discontinuities are due to finite 
size effects as explained in Sec. 12.3.11 

Transient effects and loading history. Transient flow inhomo- 
geneities may last for hours at the vicinity of the jamming tran- 
sition ll46ll : this casts doubt on band observations obtained by 
applying a given shear rate on the material initially at rest. 
Moreover, procedure-dependent steady flow inhomogeneities 
are expected at flow start-up due to normal stress differences 
trapped during the loading procedure in materials which exhibit 
significant normal stress differences lf4Sh . Much care thus has 
to be taken when analyzing inhomogeneities observed at flow 
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Figure 6: (a) Comparison between localized velocity profiles of a concentrated emulsion (squares), a foam (circles), a Carbopol gel (triangles), and a bentonite 
suspension (empty stars), (b) Dimensionless shear stress field t(R)/t v corresponding to the velocity profiles of Fig.|6jt. (c) Shear rate y vs. shear stress extracted 
from the velocity profiles of Fig. |6ji. The emulsion is composed of 1 /im diameter silicone oil droplets dispersed at a 75% volume fraction in water stabilized by 
Sodium Dodecyl Sulfate at a 8.5 wt% concentration [32]. The foam is composed of 45 /im diameter bubbles at a 92% volume fraction in a SLES foaming solution 
described in [34]. The gel is a hair gel (Vivelle Dop, France), which is mainly made up of Carbopol in water l49ll . The bentonite suspension is made up of (smectite) 
clay particles of length of order 1 /im and thickness 10 nm suspended at a 4% volume fraction in water f38Tl . 



start-up. In the Becu et al. [31] experiments, shear was ap- 
plied on a concentrated emulsion initially at rest; bands were 
observed while the shear stress was decreasing for 2 hours as 
in the Divoux et al. 14611 experiments. We did not observe a 
shear-banding behavior on the same material, when applying 
the procedure of Sec. l2.3.2l (progressive decrease of the applied 
shear rate). The bands reported by Becu et al. [31] were thus 
probably transient flow inhomogeneities. The observations of 
Rouyer et al. [41] on dry foams may also find the same ex- 
planation, as these were made when shear was applied on the 
material initially at rest, and as the experiment time was rather 
short (the applied strain was of order 1). 

Uncontrolled sources of thixotropy. A simple yield stress fluid 
can be easily turned into a thixotropic material. E.g., Ragouil- 
liaux et al. Jl5] have shown that tiny amounts of colloidal clay 
particles dispersed in the continuous phase of a simple emul- 
sion can form bridges between droplets and lead to thixotropic 
effects and to shear banding as detailed in Sec. 12.1.11 The 
same might happen with a foam. As we did not observe any 
shear banding in the same foam as Rodts et al. [33], using the 
same procedure and setup, it is suggested that the occurrence 
of shear banding in their experiments might be due to uncon- 
trolled traces of impurities in the system, e.g. clay particles, 
which were also studied with the same equipment. 



Other artefacts. Let us finally note that many observations of 
flow inhomogeneities in 2D foams are now attributed to viscous 
damping at the glass boundary $53, 54, 55 when the 

bubbles are confined by one (or two) glass plate. Wall drag 
cannot explain observations in bubble rafts with free surface; in 
these last systems, as stated above, the main possible source of 
such inhomogeneities stands in nonlocal effects. 

All results thus seem to be consistent with the absence of 
steady-state bulk shear banding in dense suspensions of soft 
particles. Of course, all of the above explanations still have to 
be checked in new experimental investigations. 

3. Consistency between macroscopic and local behavior 

In this section, we address the question of the existence of 
a single intrinsic local constitutive law to describe steady-state 
flows of yield stress fluids. If the material can be character- 
ized by such law, it is expected to account for what happens in 
any flow geometry, independently of the boundary conditions. 
Once shear banding has been dealt with, the main source of 
rheometric artefact is removed. The macroscopically measured 
behavior in rheometric experiments should then reflect the local 
behavior; it should thus a priori be trusted and used to predict 
flows in other configurations. 
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In the following, we review data from the literature where 
local and macroscopic behavior do not seem to be consistent. 
We discuss the main possible origins of this inconsistency: 
thixotropy, volume fraction inhomogeneities, and nonlocal ef- 
fects. We then present experimental results obtained in non- 
thixotropic dense suspensions of soft particles, and show that, 
under appropriate conditions, in a wide gap geometry, they are 
characterized by a single local law that matches the macro- 
scopic one. We discuss the conditions needed to obtain such 
result. 

3.1. Local observations in the literature 

The ability of the macroscopically measured stress/strain 
rate relationship to predict the flow behavior can be accurately 
tested by performing local measurements of material velocity 
profiles V(r) in geometries of controlled stress inhomogeneity 
r(r) (Couette geometry, Poiseuille flow, inclined plane flow). 
In such case, a theoretical velocity profile V,/,(r) is easily de- 
rived from the macroscopic relationship T macro = / macro (ymacro) 
as the local shear rate should be y,h(r) = /micro ( T (^); theoretical 
and experimental profiles can then be compared. An alternative 
consists in deriving the local shear rate field y(f) from the ve- 
locity profile; (r(r),y(r)) data obtained at various positions r, 
under various boundary conditions, can then be combined to 
build the locally derived constitutive behavior, which can be 
compared to macroscopic observations. 




(a) 



(b) 



Figure 7: Laws r(y) measured locally in the gap of a 1.9 cm wide gap Couette 
geometry (a) and in a w = 250 //m wide microchannel (b), during the flow of 
a concentrated emulsion of 75% droplets of 6.5 fim diameter. Different symbol 
shapes and colors correspond to different imposed rotational velocities of the 
inner cylinder (a) or different imposed pressure drops (b). The full line is a 
Herschel-Bulkley equation r = r v + ri HB y", with n = 0.5, r v = 11.6Pa and 
rj HB = 11.2Pa.s". The insets show the velocity profiles from which these laws 
are extracted. The lines in Fig. inset are fits of the velocity profiles to the 
nonlocal behavior Eq.fJ] Figure from Goyon et al. 0|. 



Such local measurements in concentrated emulsions have 
yielded surprising results |57, m some emulsions, the 

locally measured constitutive law does not match the macro- 
scopically measured one. Moreover, it seems impossible to find 



a single local constitutive law r(y) compatible with all flows. 
The apparent r(y) law needed to describe the velocity profile 
for a given boundary condition does not match those needed 
for other boundary conditions: it depends on the velocity at the 
inner cylinder in a Couette cell 11571.13 ill or the pressure gradient 
in a Poiseuille cell Jst] (see an example in Fig. I7J)). The same 
features have been observed in dense noncolloidal suspensions 
1 58 ] and 2D foams (bubble rafts) ifioll sheared in Couette ge- 
ometries. 

Three main reasons can been mentioned to explain the appar- 
ent absence of a single local constitutive law: transient effects, 
shear-induced migration of the dispersed elements, and nonlo- 
cal effects. 

3.1.1. Impact of thixotropy 

In a thixotropic material, the shear stress r is a function of 
the shear rate y and of shear history. A single relationship be- 
tween t and y can be defined only in steady state. During a 
transient, the shear rate can take many values for a given ap- 
plied stress r, depending on time and on shear history. When 
local measurements are not strictly performed in steady state, 
it may then leave the false idea that different laws are required 
to describe the flow of the sample. This is illustrated in Fig. [8] 
where local laws r(y, t) measured in a thixotropic material at 
various times t during a transient are displayed. It is clearly 
observed that the local flow behavior depends on shear history, 
until a steady-state is reached. 



ID 
CO 



30 



25 



20 



15 



, ♦ V *o ■ 

♦ ♦v v > A o°.- 

' v A O ■ 




0.1 



1 , 
Shear rate (s ) 



10 



Figure 8: Laws r(y, r) measured locally at various times t (squares: 30 s, empty 
circles: 90 s, up triangles: 150 s, empty down triangle: 210 s, diamonds: 270 s, 
full line: 1 h) during the flow of a thixotropic mud (see 1 15] for details on the 
material) in a wide gap Couette geometry. The material is first sheared at a 
rotational velocity SI = 50 rpm during 10 min; the velocity is then decreased 
abruptly at fl = 15 rpm and 2 velocity profiles are recorded per minute during 
1 h, which allows reconstructing the time-dependent local law r(y, t). 



Although dense suspensions of soft particles are usually 
found to be nonthixotropic, it has been shown that long times 
can be needed to reach a steady state at flow start-up in some 
of these systems i46ll (Sec. 12.3. U . This seems to be due to a 
specific structuration mechanism at rest; the structure formed 
at rest may then take time to be erased at low shear rate. Two 
situations should thus be distinguished: (i) when the material is 
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at rest, it gets slightly structured and its behavior at flow start- 
up may first depend on flow history during long time (as long 
as 24 h in the Divoux et al. [46] experiments); (ii) when the 
material is flowing at steady state, its subsequent flow behav- 
ior has no significant dependence on flow history as long as 
it is continuously strained in its liquid state. It contrasts with 
usual thixotropic systems such as colloidal gels, where compe- 
tition between structuration and destructuration is observed in 
the liquid state. 

These structuration effects might explain the observations of 
Becu et al. 113 ill , who cannot account for the flows of an ad- 
hesive emulsion with a single constitutive law. Their experi- 
ments are indeed performed at flow start-up. Their observation 
of flow inhomogeneity, whereas we do not observe shear band- 
ing - on the same material - when slowly decreasing the applied 
shear rate from high to low values (see previous section), then 
suggests that their local measurements are performed during a 
transient, consistently with the Divoux et al. |46f] observations. 

3.1.2. Impact of particle migration 

Particle migration can also be a cause for the apparent lack of 
a local constitutive law [17]. Indeed, when various (r(r), y(r)) 
data are combined to build a constitutive law, a hidden hypoth- 
esis is that the material is homogeneous, i.e., that the same ma- 
terial is dealt with at any position r in the flow. If there is an 
inhomogeneous volume fraction field <p{r) in the material, one 
has to be careful when inferring a constitutive law from veloc- 
ity profiles. To illustrate the problem, let us consider a volume 
fraction dependent constitutive behavior t = g(<p,y). Under two 
different boundary conditions denoted b\ and bi, a same local 
value jo of the shear rate can be reached at two positions r bl and 
r bl where the local <p values differ. When reconstructing a con- 
stitutive law from all local data, this may lead to associating two 
different stress values T b] = g((f>(r hl ),y Q ) and T bl = g((f>(r b2 ),yo) 
to the same shear rate yo, thus leading to the - erroneous - 
conclusion that flow under different boundary conditions is not 
compatible with a single constitutive behavior. 

This phenomenon explains the apparent lack of a single con- 
stitutive law to account for all flows in dense noncolloidal sus- 
pensions of rigid particles B58ll . Further studies of the same 
system have indeed shown that it is inhomogeneous [17] due 
to shear-induced particle migration. A local constitutive behav- 
ior compatible with all flows can nevertheless be built M171 |27J 
by combining (r(0(ro), y(ro)), y((/>(ro), y(ro))) data measured at 
a fixed position r*o under various boundary conditions: this pro- 
vides the constitutive behavior at a well defined volume fraction 
<p{ra) for various local shear rates y(ro). By studying this local 
behavior at various positions r*o, a ^-dependent local behavior 
is finally built. The observations of flow inhomogeneity, and 
thus of apparent viscosity inhomogeneity, in a colloidal glass 
sheared in a homogeneous stress field ll25ll is also accounted for 
by volume fraction inhomogeneities. 

This phenomenon is not expected to occur in dense suspen- 
sions of soft particles: we have shown in Sec. 12.2.21 that these 
systems remain homogeneous under shear. Consequently, it 
prob ably cannot explain the Goyon et al. ^ and Katgert et al. 
[ 10] observations. 



3.1.3. Nonlocal behavior 

The apparent absence of a single local constitutive law can 
also be the signature of nonlocal phenomena. Indeed, the local 
t(j) relationship observed in dense emulsions in microchannels 
by Goyon et al. (2D was found to depend on the boundary condi- 
tions (pressure drop, surfaces roughness), which is the hallmark 
of nonlocal behavior (Fig.[7p). 

A nonlocal modeling [9j |]JJ has been shown to successfully 
account for all flows in the Goyon et al. iS [45| experiments; 
fair agreement has also been found with the Katgert et al. ifioll 
experiments in bubble rafts. In jammed materials, flow basi- 
cally occurs locally via a succession of elastic deformation and 
irreversible plastic events. The main idea of the nonlocal model 
[59] is that a localized plastic event induces relaxation of the 
stress in the whole system, and may thus significantly affect the 
behavior in its neighborhood - its zone of influence being quan- 
tified by a length called 'flow cooperativity length'. Conse- 
quently, the rate of plastic rearrangements in a given zone de- 
pends on the rate of rearrangements in its neighborhood. It will 
then be different in a strongly inhomogeneous stress field, i.e., 
if the stress varies significantly over than in a homogeneous 
stress field far from boundaries (in this last situation, the bulk 
behavior is obtained). Nonlocal effects are thus predicted to 
take place in the presence of high stress gradients. In its sim- 
plest form, the Bocquet and Colin ill ill model is 



lyM(z) = [y/T] bM (z) + fd 2 z [j/T](z) 



(3) 



where [y/T]h u ik is the bulk apparent inverse viscosity of the ma- 
terial (also called 'fluidity') and z is the direction of the stress 
gradient |@]. This nonlocal model also accounts for the impact 
of the wall roughness, when the material display s a different 
behavior near the boundary than in the bulk 19j,|45|[. In a homo- 
geneous stress field, the length over which the behavior changes 
from surface to bulk behavior is then predicted to be of the order 
of a few £,. 

These two features explain the Goyon et al. S 45] observa- 
tions in microchannels of small size. From a fit of experimental 
data to the nonlocal model (Fig.|7p inset), they have shown that 
£ is zero below the jamming volume fraction (/>„,, and increases 
continuously with (p above (/>„,. E.g., £ is found to be 5 droplet 
diameters in a monodisperse emulsion of 85% volume fraction. 
Consistently, discrepancy between the local and macroscopic 
behavior is found in 10 to 40 droplets wide geometries but not 
in a 3000 droplets wide geometry H (see Fig. EJ. £ is found to 
be of the order of 3 bubble diameters in the bidisperse bubble 
raft studied by Katgert et al. ifioll . which is 20 bubbles wide; no 
large scale study is reported for rafts. 

In nonthixotropic dense suspensions of soft particles, the ef- 
fects discussed above can a priori be dealt with. Indeed, migra- 
tion is not an issue (Sec. I2.2.2I I and significant time-dependent 
behavior may be found only in specific conditions (Sec. 12.3. Ti l; 
moreover, the impact of nonlocal effects should be minimal in 
wide gap geometries. In the following, we propose to study 
the local behavior of such materials in well-defined conditions 
where the above mechanisms should not be at play. 
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3.2. Nonthixotropic dense suspensions of soft particles in a 
wide gap geometry 

Our goal is to see if, under some conditions, a simple lo- 
cal law accounts for yield stress fluid flows and matches their 
macroscopic response. Given the effects discussed above, ap- 
propriate setup and procedures must be used. 



3.2.1. Setup and procedure 

In order to ensure that the bulk behavior, free from the finite 
size (nonlocal) effects discussed in Sec. 13.1.31 is studied, the 
flow geometry that is used should have specific characteristics. 
The nonlocal model (Eq. [3) and the observations of Goyon et 
al. |0] first imply that the gap has to be at least tens of coop- 
erativity length £ wide, i.e., more than one hundred dispersed 
elements diameters wide to minimize the impact of the bound- 
aries. Moreover, from Eq.[3] smooth stress spatial variations (of 
a maximum order of a few % over are required. 

Here, we use the same wide gap Couette geometry as in 
Sees. 12.2.21 & 12.31 In a Couette geometry, the stress is in- 
homogeneous (Eq. Q]) and the stress relative variation over a 
length £ << Rj, where Rj is the inner cylinder radius, is 
\T(Ri + g)lT(Rd - 1| - 2£/Rj. Assuming a value of £ of the 
order of 5 dispersed elements diameter d, as found by Goyon 
et al. |@] in some systems, smooth stress spatial variations im- 
pose Rj to be at least equal to several hundreds of d. Note that 
this is only a rough estimate as these are the variations of the 
'fluidity' y/r that have to be smooth. We thus have limited our 
studies to materials whose dispersed elements diameter is less 
than 100 pm to meet the above requirements (with /?, =4.1 cm 
and gap= 1.9 cm). 

To ensure that steady-state properties are studied, the long- 
lived transient effects discussed in Sec. 12.3.11 and their pos- 
sible consequence on the local behavior (Sec. I3.1.11 l should 
be avoided. That is why we use the same procedure as in 
Sec. 12.3.21 i.e., we first impose a homogeneous flow of the 
material at high shear rate, before decreasing progressively the 
applied shear rate on the flowing - liquid-like - material (see 
Fig. Hp). 

The same materials as in Sec. 12.31 are studied: concentrated 
emulsions, Carbopol gel, and foams, described in detail in 



. 13411 . Aftera5 min preshear at 50 to 1 00 rpm (depending 
on the material), we measure both the steady-state azimuthal 
velocity profiles V{Q.,R) with MRI techniques (see Sec. 13.2. il l 
and the torque T(Q) with a Bohlin C-VOR 200 rheometer, for 
various constant rotational velocities Q of the Couette cell in- 
ner cylinder. In all systems, the stationary velocity profiles are 
found to develop within a few seconds and to remain stable 
for long durations (as long as hours in the emulsions), which 
means that these systems are not thixotropic [2]. For the Becu 
et al. 113 111 emulsion, this contrasts with the long time evolu- 
tion of the velocity profiles they observe. We argue that this 
comes from the differences in the procedures used, as discussed 
in Sec. 12.3.21 this validates the procedure we use to study a be- 
havior free from transient effects. 



3.2.2. Experimental results and analysis 

Typical steady-state velocity profiles are shown in Fig. [5] The 
constitutive laws of the materials accounting for their velocity 
profiles can be built by using both the velocity profiles and the 
torque measurements. The stress distribution t{R) at a radial 
position R within the gap is obtained from Eq.Q] we recall that 
t{Ri) = T lilnR^H) where H is the height of the inner cylinder. 
The local shear rate y{R) in the gap is inferred from the velocity 
profiles V(R) thanks to Eq. [2] Both measurements performed 
at a given radius R for a given rotational velocity Q provide a 
local data point of the constitutive law t = f(j). Note that, as 
Eqs.[T]&[2]are valid whatever the boundary conditions are, this 
analysis is not affected by a possible wall slip. As the materials 
are homogeneous, we are allowed to combine the {t(R),j(R)) 
data measured at various positions R (see Sec. 13.1.21 : they are 
a priori representative of the same material. The local constitu- 
tive laws extracted from the experiments performed at various 
rotational velocities Q on an emulsion, a foam, and a Carbopol 
gel, are plotted in Fig. [9] 
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Figure 9: Constitutive laws of (from bottom to top) a concentrated emulsion, 
a foam, and a Carbopol gel, measured locally in a wide gap Couette cell using 
MRI techniques. Macroscopic data measured in the Carbopol gel in the Couette 
geometry (crosses) and in a cone and plate geometry (filled circles) are also 
shown. The dotted lines are Herschel-Bulkley fits to the data t = t v + rj HB y" 
with: T y = 24.3 Pa, ^ Hs = 12.3 Pas 055 , and n = 0.55 (emulsion); r v = 30.3 Pa, 
T] HB = 16.2 Pas 036 , and n = 0.36 (foam); r Y = 59.3 Pa, )? H8 =29.6 Pas 042 , 
and n = 0.42 (Carbopol gel). The emulsion is composed of 6.5 /im diameter 
silicone oil droplets dispersed at a 75% volume fraction in a mixture of 50 
wt% glycerine and 50 wt% water stabilized by Brij and trimethyl tetradecyl 
ammonium bromide at a 1 wt% concentration [32]. The foam is composed of 
45 /.im diameter bubbles at a 92% volume fraction in a SLES foaming solution 
described in [34]. The gel is a hair gel (Vivelle Dop, France), which is mainly 
made up of Carbopol in water 1431 . 



A first important result is that, for each material, all the shear 
stress vs. shear rate data fall along a single curve. This means 
that, for a given material, data obtained under different bound- 
ary conditions (different Q) effectively reflect the behavior of 
a single material with a given local intrinsic constitutive law 
in simple shear. This observation contrasts with the Goyon et 
al. observations, on the same emulsion, consistently with 
the predictions of the nonlocal model discussed in Sec. 13.1.21 
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nonlocal effects are not expected in the geometry we use as the 
gap is 3000 droplets wide for this material (thus ensuring that 
both surface effects and stress inhomogeneities have no observ- 
able impact on the behavior). We also contrast our observations 
on 3D foams with those on rafts by Katgert et al. 110]; again, 
this is consistent with what can be expected from the nonlocal 
model as the gap size is 400 bubbles wide in our experiments, 
and 20 bubbles wide in the Katgert et al. Ml OH ex periments. Our 
measurements performed on the Becu et al. 113 lfl emulsion con- 
trast with their observations too. In the Becu et al. 113 lfl study, 
however, the gap being 1mm for 0.3yum droplets, the nonlocal 
effects discussed in Sec. l3.1.3l cannot be the reason for the ap- 
parent absence of a single constitutive law accounting for all 
flows. As already pointed out above, we think that the main 
difference between the two experiments stands in the procedure 
as the Becu et al. 13 ill experiments were performed at shear 
start-up on a structured material. This tends to show that, pro- 
vided an appropriate procedure is used, a single local behavior 
is recovered. 

It is worth noting that data are obtained at local shear rates as 
low as 10~ 2 s _1 , the lower measurable value depending on the 
resolution on the velocity measurement (thus on its derivative), 
and also on the coarsening rate for foams Q]. This observa- 
tion is consistent with the absence of observable shear band- 
ing (Sec. 12.31 ). All observed behaviors are finally consistent 
with a continuous local law r(y) with y — > when t — > r v . 
Good agreement is found with a Herschel-Bulkley equation 
T (t) = T y + rj HB y n for all materials in a [10~ 2 s _1 , 40 s _I ] range 
of shear rates; the fitting parameters (in particular the index n) 
depend on the material, and are given in the figure caption. 

We now get back to the macroscopic behavior. The torque 
data T(Q) obtained in the Couette geometry can be analyzed 
using standard rheometric analysis (without using the velocity 
profiles) to yield the macroscopic constitutive law r = / m aci-o(y) 
of the material; this analysis has to be performed carefully to 



account for shear localization J2 , 32 , 431 when computing the 
macroscopic shear rate from the value of the inner cylinder ro- 
tational velocity. These macroscopic data, displayed in Fig. [9] 
for the Carbopol gel, are in good agreement with the local data. 
Good agreement has also been found in the case of an emulsion 
1 32]. To push the analysis further, we have also characterized 
the Carbopol gel with a 40 mm diameter 4° angle sandblasted 
cone-and-plate geometry. In such geometry, the stress field is 
homogeneous and cannot be a source of nonlocal effects; how- 
ever, the gap still has to be wide enough to measure bulk prop- 
erties. In the used geometry, the gap varies between 1.4 mm at 
the edges and 150 //m (truncature height); the size of the dis- 
persed elements of the studied material should thus be smaller 
than « 20^m. The macroscopic data obtained in the Carbopol 
gel with this geometry are shown in Fig. [9] Good agreement 
is obtained with the macroscopic and local data measured in 
the Couette geometry. This shows that a single intrinsic local 
behavior indeed characterizes the flows of this material. This 
validates the use of macroscopic tools to study the behavior of 
these systems, provided the gap is wide enough and a procedure 
similar to the one we propose (shearing continuously from high 
to low shear rate) is used. 



Altogether, these results show that, in the strict conditions 
discussed in Sec. 13.2.11 a simple yield stress fluid behavior ex- 
ists: at steady-state, all the studied nonthixotropic dense sus- 
pensions of soft particles are indeed characterized by a single, 
intrinsic, local, continuous and monotonic constitutive law in 
simple shear. For all materials, it is moreover consistent with a 
Herschel-Bulkley equation, as usually found in standard Theo- 
logical experiments. 

4. Conclusion 

Although many complex behaviors have been observed in 
yield stress fluid flows, we have shown that, under well-defined 
conditions, a class of materials, namely nonthixotropic dense 
suspensions of soft particles, display a simple yield stress fluid 
behavior: at steady-state, they are characterized by a single, in- 
trinsic, local, continuous monotonic constitutive equation r(y) 
in simple shear with y — > when r — > t v (i.e., they are not 
shear-banding materials). 

To obtain these results, the main material characteristics are 
that: 

• they have to be nonthixotropic. More precisely, they 
should not display time-dependent behavior in their liquid 
state, otherwise they will show shear banding. Possible 
structuration at rest is not an issue if an appropriate char- 
acterization procedure is used (see below). 

• they should be made up of soft particles. Dense suspen- 
sions of soft particles have been shown to remain homo- 
geneous when sheared, in contrast to suspensions of rigid 
particles. 

Materials complying with these requirements are foams, con- 
centrated emulsions, and Carbopol gels. 

As regards the setup and procedure, two points have to be 
respected: 

• the geometry where flow is studied should have a wide 
gap, and should be characterized by smooth stress inho- 
mogeneities to ensure that nonlocal effects are not present. 
We recall that the relevant lengthscale to consider in this 
problem may be a 'cooperativity length' rather than the 
dispersed elements size, as found in concentrated emul- 
sions. 

• the following procedure should be applied: flow has first 
to be imposed at high shear rate to erase the memory of 
possible structuration at rest; afterwards, the material flow 
should be studied by progressively decreasing the applied 
shear rate, without passing by a resting period. This en- 
sures that possible transient shear banding is avoided and 
that the local flow behavior is consistent with the macro- 
scopic one. 

In these conditions, a simple yield stress fluid behavior is 
finally shown to exist. For the studied materials, the steady 
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flow behavior in simple shear can be modeled by a Herschel- 
Bulkley equation, as found in standard rheological experiments. 
A macroscopic characterization can then be trusted as it reflects 
the local intrinsic behavior of the material. It can a priori be 
used to predict the flow behavior in other configurations, pro- 
vided the studied flow also complies with the above require- 
ments (large flow scale, problem different from flow start-up at 
low shear rate). 

Of course, when confined flows or slow start-up flows are 
studied, all the complex effects discussed throughout the text 
will possibly occur; they thus still have to be studied in depth 
in fundamental studies to help predicting what happens in such 
situations. In some materials, other effects such as elastic ef- 
fects should also be incorporated in the modeling as they may 
be of crucial importance in 3D flows. 
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